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. Abstract 

in 

^ \ In supersymmetric theories with a motivated flavour structure we investigate the 

00 \ possibihty that an e' parameter as large as the measured value be generated without 

O ■ conflicting with the bounds from electron and neutron electric dipoles and /i — > 67. 

O-i. 1 Introduction 

Oh! 

D . Experiments [1] have shown beyond any reasonable doubt that direct CP violation exists in K physics, 
as conventionally described by the parameter e', expected to be non zero in the Standard Model (SM). 
A 'simple minded' comparison between published theoretical predictions for e' and experimental data 
^ . shows, however, that it is unlikely that SM physics can account for the observed value of e'. In fig. la 
we show the chronological history of the accurate theoretical predictions [2] and of the experimental 
data [1]. In fig. lb we compare the x^-probability of the theoretical and experimental values obtained 
by combining in an uncorrelated way the data shown in fig. la: the probability that a fluctuation can 
account for the discrepancy appears very small. If we believe in the data that we have employed, we 
are almost sure that new physics has been discovered. 

Even forgetting that, for various years after 1993, the experimental value of e' seemed smaller 
than the 'true' value found by the CERN experiment and confirmed by 1999 results — a fact that 
might have had some psychological influence on theorists — there is, however, at least one important 
reason for being more careful. No one of the theoretical predictions [2], except the one with the 
largest error (and the highest central value) [3] is made in a framework that includes, at the same 
time, an explanation of the 'A/ = 1/2' rule. The obvious connection between the two calculations, 
although the latter may be more difficult than the former (at least for some techniques) suggests a 
more conservative attitude with respect to the theoretical errors^. At the time being, it does not look 
implausible to think that what the experiments see is a SM effect, at least predominantly. 

This being said, it remains a relevant question if the observed e' could be produced by supersym- 
metric effects in generic and motivated supersymmetric scenarios. This question we intend to study 
in this paper. 

^For example, a value of B.R.(i3 — > Xs^) somewhat lower than the SM NLO prediction has recently stimulated a 
revision [4] of the theoretical prediction, now in perfect agreement with the experimental data. 




Figure 1: Comparison between experimental determinations (continous lines) and theoretical predic- 
tions (dotted lines) for e'/e. 



2 Supersymmetric flavour effects 

Supersymmetry can generate flavour and CP violation effects orders of magnitude above the exper- 
imental bounds (with generic soft terms) or — at the contrary — get rid of any new flavour and 
CP-violating effect (if soft terms are mediated by flavour blind interactions, like gauge interactions, 
at not too high energies). There are however two generic situations where supersymmetry gives rise 
to new effects around the experimental bounds: 

• '12/3' effects: if the sparticles of the first two generations are degenerate, while third generation 
sparticles have a different mass, 12/3 mixing angles comparable to the CKM ones are not ex- 
cluded and generate a variety of signals. Appropriately weighting possible experimental progress 
and theoretical evidence, we can try to order the signals according to their interest. Few of them 
appear rather promising 

— > 67 (together with other fi ^ e transitions), de, d^, e 
while it seems less likely (but not impossible) that effects in 

e', T fi^, Ams, b — > sss, b 57, K — > Trviy 
show up at an interesting level. 

• '1/2' effects: while generic mass matrices of Left handed and of ffight handed sfermions lead 
to unacceptably large effects (so that we assume flavour degeneracy of the first two families of 
sfermions), a generic structure of A terms gives LR mixing terms AijXij that do not unavoidably 
give too large effects. We will try to establish a hierarchy of the resulting signals also in this 
case. 

Both the 12/3 and the 1/2 effects are reasonably motivated by flavour physics: for example a U(2) 
flavour symmetry [5] can be invoked to force the degeneracy among the first two generations of 
sfermions. This symmetry does not restrict the third generation sfermions nor force the universality 
of the A terms, so that both type of effects are present in this case. 

Both type of effects may be present in unified theories with supergravity-type mediation of soft 
terms. The 12/3 effects are already there in unified theories with 'minimal' Yukawa couplings: even 
starting with universal soft terms above the unification scale, a 12/3 mass splitting between the 
sfermions is generated by the unified top Yukawa coupling [6]. The 1/2 effects require non universal A 
terms as an initial condition^ , unless one complicates the structure of the theory in a non trivial way. 

^Even small Yukawa couplings generated by operators with a more complex gauge structure than the one of the 
minimal Yukawa interactions (for example by non-renormalizable operators involving the fields that break the unified 
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The 1/2 flavour violating effects arise as follows. We assume that the Yukawa matrices for the first 
two generations contain small non-diagonal elements: for example in the d quark sector we parametrize 
them as 

dlL d2L dlL d2L 

Xs d2ij yeds 1 y ' Xs d2R \ Aueds ^22 

with similar expressions for leptons and u-quarks. The Yukawa matrix is approximately diagonalized 
by Ri2{0) rotations as Adiag = Ri2{(-sd)XR^2{^ds)- In the {d,s) basis where Adiag = diag(A(i,As) = 
diag {\edd — ^sd^ds\^ l)As is diagonal, A\ still has non diagonal complex terms 

AX = \ ^^^^'^'^dd + {M2 - A12 - A2i)esdeds] ±(^21 - ^22)es<i 

SR \ {Ai2 - ^22) fids ^22 

where ± is the sign of edd — ^sd^ds- For example, if edd = and le^^l = \eds\-, the 12 entry is proportional 
to yjuidrris and esd ~ yTTidJrris gives the dominant contribution to the Cabibbo angle. 

It has been observed in [9] that CP violating 1/2 effects generated by non universal and complex 
A-terms can give rise to a supersymmetric contribution to e' as large as the observed one. One has to 
note, however, that if the A terms were complex to start with, one would also expect a phase in the 
other soft breaking parameters (the B term, or — in the most used phase convention — the /i term) 
which would not generate e'. It is well known that electron and neutron EDMs give strong constraints 
on the phases of the ^ term, but do not disfavour a maximal phase in the A terms [10]. 

It is interesting, however, that flavour and CP violating 1/2 effects can be obtained with real soft 
terms provided that 

(1) the Yukawa couplings Xij are complex (as usually assumed in order to get a complex phase in 
the CKM matrix); 

(2) the A terms have the form AijXij with Aij real and non-universal; 

(3) there is no negligible entry in the 1,2 sector of the Yukawa matrix. 

More precisely the condition (3) means that no element of the matrix can be set to zero without 
affecting the masses in a relevant way. 

In this situation, flavour and CP-violating 1/2 effects arise in the following way. Using only 
redefinitions of the phases of the dn, dm superfields it is possible to choose a basis in which only 
one element of the 1,2 Yukawa matrix is complex. It is convenient to choose the only complex entry 
to be the 11 one {edd i^i (I))- It is now possible to obtain approximately diagonal Yukawa matrices 
by performing real rotations in the Left and i?ight 12 sectors. These rotations leave non diagonal 
elements in the ^A terms if the A terms are non universal as in eq. (2). At this point we have a 
diagonal Yukawa matrix with a complex A^ element, which has a different phase from the 11 element 
of the A\ matrix, see eq. (2). With a phase redefinition of the d^ field we can reach the 'usual' basis 
where A^ is real and positive. In this 'usual' basis the dd and the ds entries of the ^A matrix are 
complex. The same procedure can be repeated in the lepton sector. 

3 Typical hierarchy of the 1/2 signals 

We now want to study the expected hierarchy of the 1/2 signals: in particular if this scenario can be 
motivated as a supersymmetric origin of the observed value of e'. The diagonal ^-terms give rise, as 

group) are in general not enough to produce non universal A terms by RGE eflects, unlike what is claimed in [7]. We 
thank A. Pomarol and R. Rattazzi for reminding us of this point. See [8] for a discussion of decoupling effects at the 
GUT scale. 
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Figure 2: Scatter plot of values of {dj^ , e' / e) (fig. 2a) and of {de,B.R.{fi — > 67)) (fig. 2b) as generated 
by supersymmetry with non universal complex A terms and in the flavour model chosen in order to 
maximize the relative importance of e' (see text). Sampling points excluded by a too large d^ , de or 
B.R.(// — > 67) are denoted with a small dot • , while allowed points are denoted with a bigger dot •. 



usual, to electron and neutron EDMs. The non-diagonal entries give rise to e' (in the d sector), but 
also to ^ — > 67 (in the lepton sector). 

Only penguin operators are generated by 1/2 effects (so that the supersymmetric effects in e, Athb, 
K nvv' are negligible). However the hierarchy between the 'penguin effects' {fi 67, d^, d^, e' , 
and, if the analysis is extend to the third generation, b — > sss, b — > 57, r — > /X7) is not substantially 
different from the case of 12/3 effects. The main difference is that hadronic effects are not suppressed 
by possibly significant radiatively induced GIM cancellations [11]. 

We perform our analysis in a simple motivated case chosen in order to enhance the relative impor- 
tance of e' with respect to the other signals that we study. We assume that 

1. each Yukawa matrix has vanishing 11 entry and comparable 12 and 21 entries; 

2. the Aij terms are complex and non universal, while the term is real. 

With stronger simplifying assumptions on the A and AX matrices we would restrict our analysis to 
very particular situations. For example e' vanishes for antisymmetric or symmetric matrices with 
vanishing 11 element^, while the electric dipoles vanish for hermitian matrices. The zero in the 11 
position implies 1/2 mixing angles in the d-quark sector equal to eds ~ ^sd ~ v^m^/ms « 0.23, larger 
than the 1/2 mixing angles e^e ~ ^efi ~ {'me/fn^y/'^ « 0.07 between leptons. Having made these 
favourable (but reasonable) choices, we will find that a dominant effect in e' is somewhat disfavoured, 
but nevertheless possible. 

1/2 effects generate magnetic and chromo-magnetic penguins. We give in the appendix the expres- 
sions for the relevant penguin amplitudes and list here simple approximate expressions for the signals 

BR(^f ^ 67) TeV^ \A^ee^,e? + l^e^ee/.^ , . 

1.1 10-11 ^ m| 2M%me/m^ ^ ' 

^We thank Andrea Romanino and Hitoshi Murayama for having pointed out an error on this in a previous version of 
the paper 



4 



- ,.l^l^^k±J^ (3b) 



4.3 10-27ecm ' m? Mn 

dN ,TeV2 yi;^ + ^tan/3 
4 o"-*^™ — — iTT 



0.8 10-25 e cm m? M3 



676 STeV2^^e:,^:,-e,,^ 



2 10-3 2 m? 2M3y^^^^ 

In eq. (3d) the factor precisely defined in the appendix in eq. (A.l), parametrizes the uncertainty on 
the relevant iiT — > 27r matrix element. It is reasonable to expect it to lie in the range 0.3 ^-6^3. In our 
analysis we will employ the value B = 2>. A similar uncertainty affects the prediction for the neutron 
EDM, computed as described in the appendix. We have defined Xs^ds^-ds and Xs^sd-Asd as the off- 
diagonal components of the matrix in (2), and similarly for leptons. The approximate expressions 
are not valid for Mj — > 0. As discussed in the appendix, we have maybe slightly overestimated e' 
and omitted a possibly dominant contribution to the neutron EDM from the s-quark EDM [12] (this 
contribution is not present in the case of real non-universal A terms). 
The approximate expressions indicate that 

• The value of djy/e' generated by 1/2 effects does not depend strongly on the sparticle spectrum 
and is very close to the value generated by 12/3 effects (eq. 18 in [11]). With a matrix element 
for e' slightly larger than the one we have used, or with a small cancellation between the various 
contributions to c^at, the supersymmetric contribution to e' can be as large as the observed value 
without conflicting with the bound from neutron EDM. Unless for some reason the imaginary 
part of AXii in (2) is suppressed, a supersymmetric e' implies that d]\f is very close to the 
experimental bounds 

dTV = (-3 ±5)10-2^6 cm [13], c^tv = (1.9 ± 5.4)10-26ecm [14]. 



• The bound on B.R.(/i — > 67) [15] requires sleptons heavier than ~ ITeV (maybe helped by a 
light neutralino and/or by a small leptonic A term). 

• The bound from the electron EDM is not very significant. We remember that we have taken 
/i real, while the A terms have large complex phases. A complex /i term would have been a 
more efficient source of EDMs (because a large fi tan /3 is suggested in minimal models by the 
following considerations: a too low tan/3<2 tends to give a too light higgs, and a large //-term 
can compensate the large RGE gluino correction to EWSB) without enhancing e'. 

These results are confirmed by a numerical analysis, as shown in fig. 2. To study the correlation 
between hadronic and leptonic observables, we assume universal sfermion masses and universal gaugino 
masses at the unification scale. We do not require universal A terms at the GUT scale, but we assume 
that the Aij terms that multiply the d Yukawa couplings are unified with the ones that multiply the 
leptonic Yukawa couplings. We scan these parameters in the following way: 

mo, |/Uo|,M5, l^ol, 1^22!, 1^12!, \A2i\ = (0 l)msusY- 

where ttt-susy is a dimensionful parameter, and (0 1) means a random number between and 1 
(different for each parameter). The phases of the A parameters are extracted as exp[27ri(0-^ 1)], while 
the n term is taken real (positive or negative). The value of jtt-susy is determined by imposing that the 
minimum of the MSSM higgs potential be at its physical value. This scanning procedure [16] generates 
a sample of spectra with density inversely proportional to the 'fine-tuning' (i.e. configurations that 
require unlikely cancellations are more rare). 

We restrict our analysis to the few sampling spectra that still satisfy the LEP and Tevatron bounds 
on sparticle masses, that can be roughly summarized as 

mx>91GeV, mft,^85GeV, m^->80GeV M3, mg ^ (180 ^ 260) GeV. 
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Consequently in most of the cases the hghtest chargino has mass close to 100 GeV (the distribution 
of sparticle masses is plotted in fig. 4 of [16]). This procedure automatically takes into account the 
possibility of cancellations between different contributions to the observables we are studying. 

Fig.s 2 clearly show that a significant contribution to e', although possible, is not likely if we 
maintain the described correlation between quarks and leptons. This, however, may not be necessarily 
the case. 

As said, we have restricted our analysis to a simple structure of Yukawa matrix. To see what 
would happen in a different case, it is sufficient to rescale our final results by appropriate factors. For 
example it is possible to increase the value of e' /d^ without generating a too large contribution to 
the Cabibbo angle by assuming higher values of e^^. For example e^s ~ 1 can be motivated by fiavour 
symmetries. In unified models, however, this also implies e^e ~ Ij enhancing B.R.(;U — > 67) faster than 
e', see eq.s (3) (and, depending on the value of neutrino masses, giving too large Ve — > oscillations). 
This alternative choice seems thus more problematic than the one we have assumed. 

The pattern of effects generated by complex Yukawa matrices and real non universal A terms is 
qualitatively similar (although there are more unknown parameters that cannot be fixed in a simple 
motivated way), so that we do not show the corresponding plot. 

If this set of effects were observed, it will be interesting to know if they are due to 1/2 effects or 
to 12/3 effects (or to something else). Along the lines of [17], the value de/d^ could give an indirect 
discrimination of the source. If sparticles will be produced at accelerators, one can check if there is the 
12/3 non-degeneracy among sfermions (for example between and tr) characteristic of 12/3 effects. 
It might also be possible to observe fiavour violation in sfermion interactions: the amplitudes for the 
1/2 effects that we have studied will be suppressed by factors like m^/m^, while amplitudes for 12/3 
effects would only be suppressed by not very small CKM-like mixings. 



4 Conclusions 

Present theoretical uncertainties in the determination of the e' parameter do not allow a firm conclusion 
on the plausibility of explaining the measured value in the SM. In turn, the interpretation of the 
experimental result as a (main) supersymmetric effect appears highly speculative, at the time being. 

It is nevertheless true that a plausible mechanism exists [9] able to generate an effect in e'. A way 
to confirm such an effect would be to discover a neutron electric dipole just below the present bounds. 

Even with favourable structures of the Yukawa couplings, a simple minded correlation between 
quarks and leptons makes the present limit on /x — > 67 highly constraining if not excluding at all the 
possibility of a large supersymmetric contribution to e' in motivated models. This is even more true 
for the 12/3 effect than for the 1/2 effect, because of GIM-like cancellations distinguishing quarks 
and leptons in the former case. The correlation between quarks and leptons can, however, be less 
straightforward or even not exist at all. In a truly unified field theory we believe that /i — > 67 is 
anyhow the most sensitive probe of supersymmetric flavour effects. As such, the interpretation of e' as 
a (mainly) supersymmetric effect is a further motivation to look for flavour violations in the leptonic 
sector. 



A Loop effects 

1/2 effects generate magnetic and chromo-magnetic penguins. If we parametrize the coefficients of 
these operators as 



A{fL ^ f'Rl) 



;(/L7/..™^./R) + h.c. 



+ {L^R) 
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the values of the penguin amphtudes A are 



A{d^dG) = m,|5^Im^i^^^[-(2c3-C8)P;,,(r,3)-C8P;,(r,3)] 

AisL^dnG) = c8m,e:,^^[P^,(r,3)-Pij,(r,3)/9] 

A, N * eay H^^{H^^ + cot e^H^y^r) 
AinL^eRj) = -Qem^e^^^A^^^ -PsiireNj 

with analogous expressions for the u dipole and for sr di and /i/j ei transitions. The flavour 
diagonal amplitudes (/ = /') coincide with the EDMs: A{f fj) = dj and A{f fj) = dj'"^ . 
Here C3 = 4/3 and cg = 3 are QCD factors and rji = m?j/Mf is a ratio between sfermion and 
gaugino/higgsino squared masses. The P functions are defined e.g. in eq. (B.6) in [18]. If sfermions 
of opposite chirality are degenerate, P'(r) is the derivative of P{r). In general P'{rfi) stands for 
[P(rj^j) — P{t f^i)]/ {r f^i — ff^i). The neutralino mixing matrix H diagonalizes the neutralino mass 
matrix Mn as HMnH'^. 

The QCD renormalization corrections to the electric and chromo-electric dipoles have been com- 
puted integrating the RGE equations [19]^ 



dliifiUf^J 4^^ ■U?''°y' ^ V32eg,/3<73 (-29 + 2nj)/3 



from /i ~ Mz down to a QCD scale at which a^ip) ~ 47r/6. The factor raj is the number of active 
quark flavours and q = {u,d}. The neutron EDM d^ has been computed in terms of the u and d 
quark EDMs renormalized at the QCD scale using the naive quark model approximation and taking 
into account chromo magnetic dipoles using the naive dimensional analysis [20]: 

d;v-^(4d.-d.) + ^(4dQ^°-dQC°). 

We are omitting a possibly dominant contribution from the s-quark EDM [12] (not present with real 
non- universal A terms). 

The magnetic penguin operators that contribute to e' renormalize in the same way as the dipole 
operators. It is however convenient to include a g-^ in the operator (rather than in its coefficient), so 
that the anomalous dimension of the chromo magnetic dipole that gives rise to e' becomes 7 = +4/3. 
Its matrix element has been taken from [21] 

(vrvr / = 0|ff3 G^^^d\Ko) ^ ±V 2T7f ^^^T^' ~e = ^26 Re^o ^^''^ 

where, with our normalization of the meson states, F.,^ = 131 MeV, Fk = 160 MeV, Re^o = 
3.3 10~ GeV and oj = 1/22. The B factor, computed at the leading order in the chiral quark 
model, is i? PS 0.3 [21]. However B could be larger since the next orders are not suppressed by the 
factor. In our numerical analysis we have employed the value i? = 3. 
The /i ^ 67 branching ratio is given by 

B.R.(/i ^ e7) = (m3/167r)(|^(^z. ^ ejj7)l' + \A{^lR ^ eLl)\'')/T{^i). 



''We have recomputed the 7 coefhcient for dq, since its sign differs from the one in various recent phenomenological 



analyses but it agrees with the original computation [19] 
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